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I . SUMMARY 

During  the  past  several  years,  Systems,  Science  and 
Software  (S* ) personnel  have  been  actively  engaged  in  a com- 
prehensive program  involving  computer  modeling  of  the  non- 
linear processes  that  characterize  underground  nuclear  ex- 
plosions, propagation  of  the  resultant  stress  waves  through 
realistic  earth  structures  and  prediction  of  the  ground 
motion  recorded  at  teleseismic  distances  from  an  explosion. 
The  objectives  of  the  subject  project  are  to  employ  these 
modeling  and  predictive  capabilities  in  a systematic  examina- 
tion of  the  effects  of  variations  in  source  and  emplacement 
parameters  on  seismic  signals  from  underground  explosions, 
and  to  investigate  methods  for  utilizing  the  general  charac- 
teristics of  seismic  waveforms  to  obtain  reliaVjle  yield 
estimates  for  explosions. 

The  technical  phases  necessary  to  accomplish  the  ob- 
jectives of  this  project  are  as  follows; 

1.  Conduct  a systematic  theoretical  examination 
of  material,  source  and  emplacement  parameters 
which  affect  yield-magnitude  relationships 
and  compare  the  theoretical  predictions  to 
actual  observations. 

2.  Determine  and  express  uncertainties  of  yield 
estimates  in  terms  of  uncertainties  in  gross 
earth  structure,  near  source  material  proper- 
ties, and  local  source  and  receiver  structure. 

Major  accomplishments  during  the  second  three-month 
period  of  this  project  were  realized  in  several  different 
areas  of  research.  Of  particular  importance  was  the  exer- 
cise of  our  computational  capabilities  for  the  prediction 
of  teleseismically  recorded  body  and  surface  waves  from 
the  recent  Pahute  Mesa  explosion,  MAST.  Our  objective  in 
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this  experiment  was  to  model  the  near  source  (explosion) 
nonlinear  processes,  propagate  the  resultant  stress  waves 
through  realistic  earth  structures  and  finally,  generate 
synthetic  seismograms  for  comparison  with  actual  record- 
ings obtained  from  five  Special  Data  Collection  Stations 
(SDCS)  located  at  teleseismic  distances  from  the  Nevada 
Test  Site  (NTS) . 


Our  ground  motion  predictions  for  MAST  were  quite 
successful  in  terms  of  both  amplitude  and  waveform  match- 
ing. For  instance,  the  predicted  short-period  body  wave 
amplitudes  were  within  30-50  percent  of  the  observed  ampli- 
tudes at  most  of  the  SDCS  stations.  In  addition,  the 
character  of  the  first  few  seconds  of  the  P-wave  train 
at  the  various  SDCS  station  were  matched  in  some  detail. 

A technical  report  (Dache,  et  al. , 1975a)  describing  the 
initial  ground  motion  predictions  for  M7kST  was  submitted 
during  this  contractual  report  period.  A more  detailed 
report  that  describes  the  iterative  procedure  adopted  for 
the  final  ground  motion  computations  is  in  preparation  and 
will  be  submitted  shortly. 


I 


Another  major  area  of  on-going  investigation  is  the 
computation  of  the  seismic  coupling  resulting  from  the 
simultaneous  detonation  of  three  closely  "pac^d  nuclear 
sources.  Our  efforts  during  this  past  quarter  focussed  on 
the  development  of  a tension  failure  model  that  describes 
the  development  of  a region  of  enhanced  tension  failure 
(cracking)  produced  during  stress  release  behind  the  inter- 
acting shock  fronts.  This  region  is  a unique  feature  of 
the  nonlinear  material  behavior  produced  by  the  simultaneous 
detonation  of  closely  spaced  explosives.  The  multiple 
burst  calculation  has  been  carried  out  into  the  elastic 
region  using  the  two  dimensional  CRAM  code  with  the  new 
material  failure  model  ’ncluded. 
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A technical  report  (Cherry,  et  al. , 1975a)  that  des- 
cribes the  results  of  an  explosion  near  source  material 
parameter  study  was  prepared  and  submitted  during  this  re- 
porting period.  The  objective  of  the  research  presented 
in  this  report  was  to  determine  the  dependence  of  tele- 
seismic  magnitudes  on  the  nonlinear  behavior  of  the  near 
gxplosion  source  rock  environment.  Such  information  en 
2d>les  one  to  express  uncertainties  in  explosive  yield  esti- 
mates in  terms  of  uncertainties  in  the  material  properties 
and  provides  insight  concerning  the  requirements  for  col- 
lection of  geophysical  data  at  a 'specific  test  site.  The 
conduct  of  the  Limited  Yield  Test  San  Treaty  should  be 
greatly  facilitated  by  the  availability  of  this  information. 
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II.  INTRODUCTION 

As  stated  in  the  previous  section,  the  objectives  of 
the  subject  research  project  are  to  utilize  existing  compu- 
tational capabilities  to  examine  the  effects  of  various 
source  and  emplacement  parameters  on  seismic  signals  from 
\inderground  explosions,  and  to  devise  and  evaluate  methods 
for  utilizing  the  general  characteristics  of  seismic  wave- 
forms to  obtain  reliable  estimates  of  explosion  yields.  In 
order  to  realize  these  objectives,  activity  on  this  program 
during  the  second  three-month  period  of  this  contract  con- 
centrated in  the  following  areas: 

1.  The  prediction  of  teleseismic  ground  motion, 
body  emd  surface  waves,  generated  by  specific 
NTS  explosions  and  recorded  at  selected  seismo- 
graph stations. 

2.  Modeling  of  the  source  region  of  a multiple 
underground  nucleau:  explosion  scenario  and 
computation  of  the  seismic  coupling  from 
such  an  event. 

3.  Investigation  of  the  dependence  of  the  seismic 
coupling  of  a nuclear  explosion  on  the  non- 
linear behavior  of  the  near  explosion  source 
rock  environment. 

The  plan  of  the  remainder  cf  tJiis  report  is  to  pre- 
sent technical  discussions  of  each  of  these  three  research 
areas,  followed  by  a section  summarizing  the  most  important 
results  obtained  to  date. 
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III.  TECHNICAL  DISCUSSION 

3.1  GROUND  MOTION  PREDICTIONS  FOR  THE  MAST  EVENT 

It  was  requested  that  S*  predict  the  short  and  long 
period  seismograms  recorded  for  the  Mast  event  at  five  specified 
SDCS  stations.  An  initial  prediction  was  made  using  the  best 
available  information  for  the  rock  properties  at  the  working 
point  and  the  near  source  geology.  Standard  earth  models  were 
used  for  the  propagation  path.  Upon  comparing  our  predictions 
to  the  actual  data,  we  concluded  that  much,  if  not  most,  of  the 
observed  discrepancy  was  due  to  propagation  path  effects. 
Improved  modeling  of  the  propagation  path  is  expected  to  lead 
to  improved  predictions  of  future  events. 

In  this  section  we  first  summarize  our  initial  predic- 
tions of  Mast  and  the  comparison  with  data.  After  studying 
the  comparison  we  made  adjustments  to  our  models  of  the  propa- 
gation path  and  the  improved  results  are  also  summarized. 

3.1.1  Body  Waves 

The  computational  method  for  the  prediction  of  body 
waves  employs  the  following; 

1.  The  reduc'jd  displacement  potential  (RDP)  is 
computed  by  our  spherically  symmetric,  one- 
dimensional  finite  difference  code  (SKIPPER) . 

2.  The  RDP  is  input  to  a code  which  computes  the 
detailed  crustal  reverberations  for  an  equiva- 
lent elastic  source  buried  in  a stack  of  plane 
elastic  layers  (see  Appendix  B,  Bache,  et  al. 

(1975b)  for  the  theory) . 

3.  Detailed  crustal  reverberations  at  the  receiver 
are  computed  using  the  method  of  Haskell  (1962) . 
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4.  Travel  through  the  upper  mantle  is  computed 

t by  generalized  ray  theory  aa  implemented  by 

Helmberger  and  described  in  V^iggins  and 
He  Lmberqer  (1974). 

5.  The  desired  sensing  instrument  transfer 

^ function  is  applied  to  the  ground  motion. 

For  surface  waveS/  steps  2-4  are  replaced  by: 

2.  The  Rayleigh  wave  dispersion  data  and  amplitude 

* excitation  for  a plane  layered  earth  model  is 

computed  using  the  method  of  Harkrider  (1964/ 
1970) . 

3.  The  source  time  history  is  convolved  with  the 
earth  response  function. 

The  pertinent  data  were  obtained  from  the  fol]  owing 
sources: 

1.  The  material  properties  of  the  source  region 
were  extracted  from  the  CEP  data  synopsis  of 
6/18/75  (Cherry,  et  al./  1975b). 

2.  The  velocity-depth  profile  for  the  region  above 
the  working  point  was  deduced  from  the  same 
source.  We  also  obtained  a log  to  a depth  of 
2.3  km  of  the  nearby  hole  UE19d  from  Dick  Ramspot 

f of  LLL.  This  was  used  to  construct  a profile 

to  that  depth.  The  remainder  of  the  crustal 
model  for  the  body  wave  calculations  is  based 
on  an  average  basin  and  range  crust  of  Hill 
^ and  PakJ.ser  (1967) . 

3.  since  we  had  little  or  no  information  about 
crustal  structure  at  each  of  the  receiver 
locations  specified  by  the  Project  Officer, 

^ the  Scune  average  crustal  model  was  used  at  all 

stations  for  the  body  wave  calculations. 


8 


R-2788 


4.  For  the  initial  prediction  of  body  waves,  two 
upper  mantle  models  were  used.  These  are  KWNE 
from  Helmberger  and  Wiggins  (1971)  and  an 
unpublished  model  by  Anderson,  Hart  and  Jordan 
called  C2AJ.  Detailed  studies  of  data  re- 
corded at  stations  on  a northerly  azimuth  and 
30-36*  from  NTS  found  C2AJ,  slightly  modified, 
to  be  the  preferred  model.  At  the  time  of  the 
first  prediction  we  had  not  accumulated  enough 
experience  at  other  azimuths  and  distances  to 
make  a reliable  selection  of  an  appropriate 
model. 

5.  For  the  initial  prediction  of  surface  waves 
the  earth  model  CIT109  (Archambeau,  Flinn  and 
Lambert  (1969))  was  used,  mainly  because  it 
was  computationally  convenient.  However,  this 
model  was  primarily  based  on  inversion  of  body 
wave  data  and  the  crustal  layers  were  meant  to 
be  generally  applicable  to  shield  areas. 

6.  For  the  revised  predictions  of  body  waves,  all 
features  of  the  calculation  were  held  the  same 
except  different  earth  upper  mantle  models  were 
used.  The  new  models  were  a modified  version 
of  HWNE  and  a modified  version  of  HWA  (Wiggins 
and  He-mberger  (1973)).  The  modifications  were 
primarily  motivated  by  a study  of  the  Mast 
observations. 

7.  We  are  currently  evaluating  a number  of  crustal 
models  for  an  improved  prediction  of  surface 
waves. 

8.  The  instrument  transfer  functions  were  those 
provided  by  the  Project  Officer. 
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Incidental  to  the  predictions,  we  have  been  evaluating 
the  quality  Uor  short  period  observations)  of  the  SDCS  sta- 
tions at  which  the  predictions  are  being  niade.  The  station 
quality  is  quantized  in  t#^x.is  of  (1)  the  amount  of  energy 
converted  to  tangential  ground  motion;  (2)  the  coda  length 
which  is  representative  of  internal  refleouion  and  scattering 
at  the  receiver,  and  (3)  the  amplitude  coupling  of  the  ptation. 
The  latter  is  a quantization  of  that  poir.ion  of  the  amplitude 
response  of  a station  which  is  independent  of  azimuth  and 
distance. 

The  first  conclusion  from  our  station  quality  estimate 
is  that  the  FNWV  site  is  of  low  quality  and  its  usefulness 
(for  short  period  observations)  to  the  SDCS  netv.or)c  should 
be  reevaluated.  We  need  to  obtain  more  data  to  firmly 
evaluate  the  other  stations  but  are  in  the  process  of  doing 

so. 

The  comparison  between  our  theoretical  and  the  observed 
seismogrcims  will  now  be  presented.  In  Figs.  3.1  and  3.2  the 
initial  short  period  predictions  are  shown  together  with  the 
observations.  Computed  seismograms  for  the  two  upper  mantle 
models  C2AJ  and  HWNE  are  shown.  From  both  the  theoretical 
find  observed  records  two  amplitudes  were  measured.  The 
first  is  the  "b"  amplitude  which  is  the  peak-to-peak  amplitude 
of  the  first  cycle  on  the  record,  corrected  for  the  instrument 
response  at  the  apparent  period  (i.e.,  twice  the  time  separa- 
tion between  the  measured  peaks) . The  second  is  denoted  the 
"d"  phase.  The  nomenclature  arises  from  the  convention  for 
measuring  m^^  followed  for  the  LRSM  network  herein  the  maximum 
cycle  in  the  first  few  cycles  on  the  wave  train  is  measured. 
The  cycle  used  for  this  amplitude  is  indicated  on  each  seismo- 

grcim. 

The  comparison  between  theoretical  and  observed  "b" 
and  "d"  phases  from  Figs.  3.1  and  3.2  are  shown  in  Figs.  3.3 
and  3.4.  Excimining  the  wave  form  comparison  in  Figs.  3.1  and 
3.2  we  conclude  that  many  of  the  important  features  of  the 
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seismogram  are  predicted  by  the  model  but  that  th>  compari- 
son is  far  from  ideal.  The  theoretical  b eunplitudes  are  sub- 
stantially smaller  than  the  observed.  Our  models  should  be 
most  accurate  in  predicting  this  portion  of  the  wave  train. 

The  agreement  between  d amplitudes,  the  phase  from  which 
conventional  m^^  is  determined,  is  somewhat  better.  However, 
this  is  to  some  extent  fortuitous  since  all  the  interfering 
phases  which  contribute  to  the  d amplitude  do  not  appear  to  be 
properly  included  in  our  calculations. 

An  effort  was  made  to  select  2m  Improved  upper  mantle 
model  to  give  closer  agreement  between  theory  and  observa- 
tions. Implicit  in  this  exercise  is  the  assumption  that  the 
source  2md  its  vicinity  is  correctly  modeled.  A separate 
study  (Bache,  ^ ^1. , 1975c)  shows  that  such  confidence  is 
warranted. 

The  results  of  the  more  recent  theoretical  seismo- 
gram calculations  are  shown  in  Figs.  3.5  and  3.6  for  two  earth 
models:  HWA-2  and  HWNE-3.  The  amplitude  comparisons  for  the 

"b"  and  "d"  amplitudes  are  shown  in  Figs.  3.7  and  3.8.  In 
our  study  of  station  characteristics,  the  FNWV  station  was 
found  to  give  amplitudes  that  are  three  tines  smaller  than 
expected,  irrespective  of  epicentral  distance.  Therefore, 
the  observed  2unplitudes  for  FNWV  plotted  in  Figs.  3.7  and 
3.8  have  been  multiplied  by  the  station  correction  factor  of 
three. 

The  waveform  comparisons  for  the  new  predictions. 

Figs.  3.4  and  3.5,  are  perhaps  a bit  better  than  for  the 
.nitial  predictions.  Figs.  3.1  £md  3.2.  Our  model  still 
does  not  successfully  account  for  all  the  early  arriving, 
interfering  phases  that  control  the  shape  of  the  first  few 
seconds  of  the  RKON  and  CPSO  records.  However,  it  does 
appear  that  the  factors  controlling  the  amplitude  of  the  b 
phase  are  properly  included  in  the  model.  Evaluation  of  the 
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predictions  should  rely  most  heavily  on  the  b phase  compari- 
son since  this  phase  is  least  affected  by  phenomena  not 
related  to  the  explosion  coupling  into  elastic  waves. 

For  the  b 2unplitude  the  comparison  (Fig.  3.6)  be- 
tween predicted  and  observed  data  is  quite  favorable / parti- 
cularly for  the  HNME-3  model.  At  several  stations  the  two 
models  span  the  observations.  Taking  HNME-3  aloi.e,  the 
eunplitudes  are  consistently  within  30-50  percent  c I the 
observations. 

For  the  d amplitude,  the  measurements  are  much  more 
scattered,  as  might  be  expected. 

3.1.2  Rayleigh  Waves 

The  comparison  between  predicted  and  observed  Rayleigh 
surface  waves  for  the  propagation  path  model  CIT109  is  shown 
in  Figs.  3.9  and  3.10.  The  amplitude  measurements  in  Fig.  3.10 
are  Airy  phase  measurements  made  on  the  indicated  cycle  on  the 
long  period  vertical  (LPZ)  seismogr2uns.  The  amplitudes  are 
corrected  for  instrument  response  at  the  apparent  period  of 
the  measured  cycle. 

As  is  immediately  ?pparent  from  the  observed  and  pre- 
dicted travel  times,  the  model  CIT109  is  not  very  appropriate 
for  these  stations.  Since  the  predicted  travel  times  are  much 
lower  than  the  observed,  the  crustal  velocities  in  the  model 
must  be  too  high.  As  far  as  the  ariplitude  comparison  in 
Fig.  3.10  is  concerned,  the  predicted  Airy  phase  amplitude 
agrees  quite  well  with  the  observed  at  RKON,  a station  on  the 
Canadian  shield  for  which  the  model  CIT109  is  perhaps  most 
appropriate.  At  other  stations,  however,  the  predicted 
eunplitudes  are  too  low  by  an  average  factor  of  =2.  We  are 
currently  constructing  more  appropriate  propagation  path 
models  and  will  construct  revised  theoretical  records  which 
should  exhibit  much  better  agreement  with  the  data. 
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3.2  MULTIPLE  BURST  CALCULATION 

3.2.1  Introduction 

Tension  failure  models  employed  in  the  CRAM  code  have 
provided  insufficient  information  for  many  applications  be- 
cause crack  strains  were  calculated  entirely  on  the  basis  of 
the  current  magnitude  and  orientation  of  the  principle  co- 
ordinates. The  crack  direction  was  allowed  to  rotate  with 
the  principle  coordinate  system  and  therefore  no  meaningful 
information  about  crack  orientation  was  provided.  A new  ten 
sion  failure  model  which  allows  for  the  adjustment  of  crack 
strain  components  in  a system  of  "failure  coordinates"  while 
holding  the  principle  direction  invariant  is  described  in 
this  section.  This  new  model  has  been  tested  in  a jCRAM  cal- 
culation involving  simultaneous  detonation  of  three  closely 
spaced  underground  explosions  and  a graphical  display  j.s 
used  to  vividly  depict  the  effect  of  shock  interaction  on 
crack  orientation. 

3.2.2  Discussion  of  the  Tension  Failure  Model 

An  element  is  assumed  to  fail  in  tension  if  a principle 
stress  is  greater  than  zero  and  if  the  element  has  at  any  time 
experienced  shear  failure.  We  then  apply  the  tension  failure 
model  proposed  by  Maenchen  and  Sack  (1964)  and  introduce  an 
inelastic  strain  normal  to  the  crack.  This  inelastic  strain 
is  just  sufficient  to  zero  the  principle  tensile  stress.  The 
magnitude  of  the  crack  strains  may  be  adjusted  at  later  times 
to  account  for  changes  in  the  stress  field,  but  the  initial 
orientation  of  the  crack  is  changed  only  to  account  for  actual 
rotation  of  the  element. 

For  excimple,  in  a cell  which  has  experienced  shear 
failure  but  remains  uncracked,  if  a , a and  o are  the 

^11  22  33 

three  principle  stresses  and  if  a is  greater  than  zero. 
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then  the  adjusted  stresses  (a  , a , a ) are  given  by 

112211 


**a  + 

11  11  \ J / 11 

a -a  -/k-^vilAe 

22  22  \ ^ I 

*=a  -/k-x-vi^Ae 
11  11  \ 3 / 11 


where 


(3.1) 


(3.2) 


k + f u 

k is  the  bulk  modulus,  vi  is  the  shear  modulus  and  all 
stresses  include  the  overburden  pressure. 

If  two  of  the  principle  stresses,  say  a and  o , 

11  2 2 

are  greater  than  zero,  then  the  stress  adjustment  becomes 

a «a  -(k-|*’-|  (Ae  + Ae  ) - 2u  Ae  , 

11  1 1 \ ^ ^ • 22  11 

a =0  “ (k  - 4 (Ae  Ae  ) - 2y  Ae  , (3.3) 

22  22  \ 3/  11  22  22 

a =0  -(k-|-y)  (Ae  + Ae  ) , 

31  33  \ 3/  11  22 


a =0 

3 3 3 
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The  above  procedure  is  followed  only  for  the  initial 
formation  of  the  crack.  The  orientation  of  the  principle 
coordinates  on  the  cycla  during  which  a crack  forms  is  de~ 
fined  by  the  angle  a (see  Fig.  3.11)  which  is  adjusted  only 
for  actual  rotation  of  thf,  element  during  the  remainder  of  a 
calculation.  The  coordinates  located  by  the  angle  a are 
defined  as  the  "failure  coordinates"  (I",  2',  2').  All  sub- 
sequent adjustments  to  the  crack  strains  are  performed  in 
the  failure  coordinate  system.  The  tensile  crack  strains 
are  adjusted  such  that  all  tensions  which  appear  in  the  element 
are  eliminated.  Additionally,  a shear  component  of  the  crack 
strain  v.ensor  is  introduced  and  adjusted  such  that  the  orien- 
tation of  the  principle  coordinates,  (1,  2,  3),  remains 
constant.  The  relative  orientation  of  the  (X,  Y) , (2,  3)  and 
(2",  3'*)  coordinates  is  illustrated  in  Fig.  3.12.  Note  that 
the  principle  coordinates  are  located  by  the  angle  (j)  and 
that  0 **  ci-<^  is  the  angles  between  the  principle  coordinates 
and  the  crack  coordinates. 

For  example,  in  a cell  which  has  already  failed  in 
tension,  if  and  are  the  failure 

stresses  and  if  o ^ ^ is  greater  than  zero,  then  the  adjusted 
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Figure  3.12.  Relative  orientatior  of  the  principal  and  failure  coordinate  systems 
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where 
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k + T U 
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Ae  = 
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- Ae 


JJli. 


tan  29  “ 0 


If  the  failure  stresnes  a and  a are  both 

l'l“  2''2" 

greater  than  zero  then  the  stress  adjustment  becomes 


2y  Ae 


-fk-i-Mj  (Ae  +Ae  )- 

1*^1'  l"!"  \ 3 / l'l“  2 "2" 

(Ae  + Ae  ) - 2y  Ae 
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All  the  normal  inelastic  strain  increments  (Ae 
Ae  ^ Ae  ^ are  accumulated  on  each  cycle  giving 


• 4 


R-2788 


„n+l  , j.n  + , 

I'l'  I'l'  I i 

n+1  « e"  + Ae  , 

^2'!'  2'2'  2'2' 

3't'  3'3' 

These  equations  apply  equally  well  for  crack  closure.  If  at 
least  one  of  the  total  strains  in  Eq.  (3.9)  is  greater  than 
zero  (say  e"  ) , then  the  crack  will  open  or  close  dependi.ng 

on  the  sign'of  If  ^ ° **'®" 

Ae  > 0 » 

1 1 

n+1  > e’^^  ^ , 

I'l'  11 

and  the  crack  width  increases.  The  inequalities  are  reversed 

df  5 <0  and  the  crack  width  decreases.  Closure  will 

1 "i 

continue  until 


e’''  + Ae  ^ < 0 • 

I'l'  I'l' 


Then 


E"t^  * 0 . 

I'l 


and  the  crack  is  completely  closed.  When  this  state  is 
achieved  the  element  is  then  able  to  support  a compressive 
Stress  in  the  (l*^#  0/  0)  direction. 
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3.2.3  Results  of  Multiple  Burst  Calculation 


The  new  tension  failure  model  has  been  incorporated 
into  the  CRAM  code  and  used  in  the  calculation  of  ground 
motion  due  to  simultaneous  detonation  of  three  15  kton  ex- 
plosive sources.  The  sources  were  spaced  165  meters  apart 
on  a horizontal  axis.  The  calculation  was  initiated  with 
the  SKIPPER  code  in  spherical  symmetry  and  linked  to  the 
CRAM  in  axial  symmetry  by  initializing  the  CRAM  mesh  just 
prior  to  shock  interaction  between  adjacent  sources. 


Grid  plots  showing  the  formation  of  cracks  as  shock 
interaction  occurs  are  presented  at  successive  times  in 
Figs.  3.13,  3.14,  3.15  and  3.16.  The  right  boundary  in 
these  figures  is  a plane  of  symmetry  and  the  lower  boundary 
is  an  axis  of  symmetry.  Elements  which  have  developed  cracks 
are  marked  with  dots  and/or  lines.  The  dots  represent  cracks 
which  lie  in  planes  passing  through  the  axis  oi:  symmetry. 

The  lines  represent  cracks  which  lie  in  planes  which  are 
perpendicular  to  the  figure  and  pass  through  the  given  line; 
thus  the  lines  provide  a convenient  display  of  the  crack 
orientation.  Note  that  in  Fig.  3.13  hoop  cracks  due  to 
radial  expansion  of  the  cavities  are  beginning  to  form. 

In  Fig.  3.14  shock  interaction  has  caused  adjacent  sources 
to  be  connected  by  a series  of  hoop  cracks.  Extensive 
cracking  appears  above  and  between  the  adjacent  sources  in 
Fig.  3.15;  this  cracking  is  in  the  plane  of  the  grid  plot 
and  indicates  that  the  two  spherical  sources  resemble  a line 
source  at  moderate  distances  from  the  axis  of  symmetry. 

In  Fig.  3.16  cracks  have  covered  most  of  the  grid. 


3.2.4  Summan 


A new  tension  failure  model  has  been  developed  and  in- 
corporated into  the  CRAM  code.  The  new  model  retains  the 
orientation  of  the  crack  at  initial  failure.  Subsequent 
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adjustment  o;;  the  crack  strains  is  accomplished  in  a manner 
which  holds  the  principle  direct ‘.ons  invariant  thereby 
preserving  stability. 

3 . 3 SEISMIC  COUPLING  FROM  A NUCLEAR  EXPLOSION 

3.3.1  Dependence  of  Seismic  Coupling  on  the  Near  Explosion 

Source  Environment 

A special  report  that  describes  a parameter  study  of 
the  dependence  of  seismic  coupling  on  the  nonlinear  behavior 
of  the  near  explosion  source  environment  was  submitted  during 
this  reporting  period  (Cherry/  e^  al. / 1975a) . The  results 
of  this  parauneter  study  provide  valuable  insight  concerning 
the  importance  of  various  material  parameters  with  respect 
to  accurate  prediction  of  seismic  coupling.  In  this  study 
it  was  shown  that  increasing  the  overburden  pressure,  or, 
equivalently,  the  depth  of  burial,  decreases  coupling  effi- 
ciency substantially.  Therefore,  the  DOB  must  be  accurately 
known  in  order  to  makr*  a reliable  prediction  of  seismic 
magnitude.  It  is  also  desirable  to  know  the  properties  of 
the  material  in  the  immediate  vicinity  of  the  device;  if 
substantial  spatial  variation  of  the  material  properties 
occurs  near  the  device  or  if  the  DOB  is  quite  shallow  it 
may  be  necessary  to  consider  two-dimensional  effects  which 
were  not  treated  in  the  pa. imeter  study.  Accurate  predic- 
tion of  seismic  magnitudes  via  one -dimensional  spherically 
symmetric  calculations  is  therefore  understood  to  be 
dependent  on  the  material  uniformity  about  the  device. 

In  this  study  a relationship  between  teleseismic 
magnitude,  rrij^,  and  '?(«),  the  steady  state  value  of  the 
reduced  displacement  potential  (RDP) , was  given  as 


itij^  log  [c  y {<”)  ] 
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where  c is  the  near  source  compressional.  wave  speed. 
Based  on  Eq.  (3.10),  if  yi(“)  and  are  two  values  of 

y (oo)  corresponding  to  materials  i and  k then  the 
changes  in  teleseismic  magnitude  is  given  by 


Am  “ 


log 


Ci 


(3.11) 


Since  the  elastic  properties  of  the  near  source 
material  appear  in  the  magnitude  relation,  Eq.  (3.11),  via 
the  P wave  velocity,  c,  it  is  obviously  important  to  deter- 
mine their  values  accurately.  A positive  error  in  deteirmina- 
tion  of  c would  cause  body  wave  magnitude  to  be  over- 
predicted. A positive  error  in  u would  cause  magnitude 
to  be  under-predicted. 

The  steady  state  value  of  the  RDP,  '!’(»),  which  appears 
in  the  magnitude  relation,  Eq.  (3.11),  is  dependent  on  the 
shock  response  of  the  near  source  material.  The  near  source 
material  properties  which  have  the  most  pronounced  effect 
on  the  shock  response  are  the  water  mass  fraction,  f^,,  and 
the  air-filled  void  fraction,  Positive  errors  in  either 

f or  (|)  could  lead  to  substantial  under-prediction  of 
seismic  magnitudes.  Seismic  magnitude  is  not  very  sensitive 
to  P , the  elastic  pressure,  and  P^,  the.  crush  pressure,  indi- 
cating thSt  details  in  the  porosity  model  are  relatively 
unimportant. 

Seismic  magnitudes  are,  however,  very  sensitive  to 
the  paraumeters  which  describe  the  failure  surface.  If  Y^, 

Y , or  P are  varied  such  that  the  material  strength  is  en- 
hcmced  the  coupling  efficiency  of  an  explosive  device  is  im- 
paired. Thus,  a positive  error  in  the  material  strength 
Igcid  one  to  under —predict  seismxc  magnxtude. 
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3.3.2  Constitutive  Modelin< 


A topical  report  by  Garg  (1975)  was  prepared  and 
submitted  during  this  contractual  report  period.  This  re- 
port deals  with  the  development  of  constitutive  re3  ations 
for  fluid-saturated  porous  media,  suitable  for  incxusion 
in  the  CRAM  or  SKIPPER  codes.  The  theoretical  formulation 
is  based  on  the  models  for  fluid-saturated  rock  aggregates 
previously  developed  by  Garg  and  Nur  (1973) . A method  for 
including  the  new  constitutive  model  in  the  standard  hydro- 
dynamic  codes  referred  to  above  is  also  described  in  this 
report. 
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IV.  CONCLUSIONS 


Considerable  progress  was  made  during  the  second 
three-month  period  of  this  project  on  several  of  the  tasks 
specified  in  the  work  statement  of  this  contract. 

Teleseismic  ground  motion  predictions  for  the  recent 
Pahute  Mesa  explosion,  MAST,  were  quite  successful  in  terms 
of  both  amplitude  and  waveform  matching.  The  predicted 
short-period  body  wave  amplitudes  were  within  30  to  50  per- 
cent of  the  observed  aunplitudes  at  most  of  the  SDCS  stations. 
In  addition,  the  general  character  of  the  first  few  seconds 
of  the  P-wave  trains  at  the  various  SDCS  stations  were 
matched  in  reasonable  detail. 

A tension  failure  model  that  describes  the  develop- 
ment of  a region  of  enhanced  tension  failure  (cracking) 
produced  during  stress  release  behind  the  interacting 
shock  fronts  from  closely  spaced  explosions  was  developed. 
Calculations  for  a multiple  explosion  scenario  (three 
closely  spaced  explosions  detonated  simultaneously)  were 
carried  out  into  the  elastic  region  using  the  two-dimen- 
sional CRAM  code  with  the  new  material  failure  model  in- 
cluded. 
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